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Abstract. We present the results of exploratory HST 
observations to detect the emission lines of ions in the 
P Pictoris disk with the spectrograph slit placed off the 
star. Possible emission lines from Fe ii have been detected 
at ~ 10~^'*erg cm~^ s~^ at 0.5" from /? Pictoris, which 
would suggest a total Fe ii density of ^ 2 • 10~^cm~^ at 
10 AU. This detection needs confirmation. If real, it re- 
quires a large production rate of gas and dust equivalent 
to the total disruption of ten bodies 30 kilometers in ra- 
dius per year. 

Key words: stars: (3 Pic - circumstellar matter - plane- 
tary systems 



1. Introduction 

The /3 Pictoris disk is an evolved replenished circumstellar 
disk around a main sequence star, and appears to be a 
planetary (or cometary) disk (Vidal-Madjar et al. 1998). 
In addition to the dust disk, gas has been probed through 
absorption lines on the star line of sight (Lagrange et al. 
1998). The observations of this gas component histori- 
cally gave the most unexpected and important results with 
the detection of what appears to be the first extra-solar 
comets ever observed (Ferlet et al. 1987, Vidal-Madjar et 
al. 1998). However, non-radial motions are still unknown 
because observations are limited to the gaseous absorp- 
tion lines on the star line of sight. Here we present the 
results of exploratory HST observations aimed to observe 
the emission lines of ions in the disk with the spectrograph 
slit off the star. 



2.1. Observations 

The observations were made on February 9, 1996 with the 
HST/GHRS using the SSA slit (0.22"x0.22") and echelle 
spectroscopy at an expected resolution of 80 000 to get the 
highest chance of detection of lines with a priori unknown 
width. After a peak-up on the star /3 Pictoris and the ac- 
quisition of reference spectra, the telescope was shifted to 
the south-west part of the /3 Pictoris disk at a position an- 
gle of 31.5 degrees to a distance of 0.5 and 1.5" from the 
star. The log of the observations is summarized in Table |[ 
The candidate lines are those from ions for which strong 
absorptions have been seen in the stable component due 
to the gas disk and in the variable components due to the 
infalling comets. On this basis, we expected the strongest 
fluorescence. 



Table 1. Log of the observations performed on Feb. 9' 
1996. r is the distance from (3 Pictoris 



# 


r 


Wavelength 


Lines 


Time 


Texp 






(A) 




(UT) 


(s) 


1 


0.0 " 


1803.1 - 1812.7 


Ni II 


08:39:03 


230 








S I, Si II 






2 


0.0 " 


2373.0 - 2384.6 


Fe II 


08:46:42 


173 


3 


0.5 " 


2373.0 - 2384.6 


Fe II 


08:55:45 


2016 


4 


0.5 " 


1803.1 - 1812.7 


Ni II 


10:28:45 


2073 








S I, Si II 






5 


1.5 " 


1803.1 - 1812.7 


Ni II 


11:58:39 


3571 








S I, Si II 






6 


0.5 " 


1652.7 - 1661.2 


C I 


14:46:39 


2534 








C I*, C I" 







2.2. The method 



2. Observations and data analysis 



The spectra observed on the disk present the same gen- 
eral shape as the spectra on the star but with a lower 
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intensity level. Indeed, any spectrum of the disk obtained 
off tlie star (Fobs (A)) is an addition of a "noise" compo- 
nent and a spectrum carrying information from the disk 
at the pointed position. The "noise" is the spectrum of 
the starhght scattered by the telescope (ai^*(A), where 
Ff{X) is the star spectrum) plus a background level (6) 
due to the fact that the background determined by the 
HST pipeline can be slightly miscalculated. The part of 
the spectrum due to the disk {Fdisk) is the starlight scat- 
tered by the dust and the possible emission lines due to 
gaseous fluorescence. We thus consider that the additional 
noisy component {Fnoise) is a linear combination of the 
stellar spectra (F*(A)) obtained at the same wavelength a 
few orbits before: 

Fobs{\) = Fno^se{\) + Fdisk{\) ~ aF^{\) + b + FdiskW{l) 

To put forward the potential presence of a component dif- 
ferent from the dominant starlight scattered by the tele- 
scope, we divide the observed spectra by the stellar spec- 
tra. The presence of a disk contribution at a wavelength 
A will be detected if the ratio Fofcs(A)/(ai^,(A) -I- b) is sta- 
tistically significantly different from 1. 

The main problem is thus the determination of a and 
b, and their errorbars within a confidence level. Assuming 
that the disk spectra (dust scattered light plus gas emis- 
sion lines) have a negligible contribution to the observed 
spectra (Fnoise (A) ^ FdiskW), wc can determine a and b 
by a minimization of 

= $]wA.(aF,(A,) + &-Fofcs(A,))', (2) 

where w\. — ^/<j\. is the weight of each measurement at 
Aj . This procedure gives not only the best determination of 
a and b but also intervals of confidence for these constants. 

2.3. Results in the Fe ii lines 

Spectra in the Fe ii wavelength range were obtained less 
than 10 minutes apart with the same instrument setting. 
The star and the disk spectra can be easily superimposed 
as shown in Fig. |l|. An excess of flux clearly appears in 
the blue and the red part of the two strongest Fe II lines 
at rest wavelength 2374.461 A and 2382.765 A. If a and 
6 are determined as described in Sect |2.2| , we get = 
197 ± 3 and 6 = (3 ± 2) • lO^^^ erg cm^s^^ A^^. With 
these parameters, a plot of the ratio of the disk to the 
star spectra reveals an excess emission at 2a level (Fig. 
The flux from the disk can be evaluated by the difference 
between the two spectra (Fo6s(A) — (aF*(A) -I- 6)), and 
is Fdisk ~ 10""'^^ erg cm~^ A~^ between 50 and 150 
kms^^ in the red part of the lines and around —50 kms^^ 
in the blue part (Fig. ||). This emission is about 1 A wide 
and stronger in the red part of the lines. 



2.5x10-13 m — ' — ' — ' — ' — I — ' — ' — ' — ' — \ — ' — ' — ' — ' — I 5x10 11 
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Fig. 1. Plot of the disk spectrum (thin line, left axis) and 
the star spectrum (thick line, right axis) in the region of 
the Fe ii hues at 2374.461 A (top panel) and 2382.765 A 
(bottom panel) . The velocities are relative to (3 Pictoris 
(21 kms^^ heliocentric). The data on the disk have been 
rebinncd by 5 pixels, the data on the star by 2 pixels. The 
ratio of the right axis to the left axis is 197. 

2. 4- Discussion 

The detection of apparent excess emission in the two Fe ii 
lines can be explained in different ways. We propose that 
this can be a real detection of the emission through the 
scattering of the starlight by the Fe ii ions in very high ve- 
locity motions. This will be discussed in detail in the next 
section. However other possibilities must be evaluated. 

First, this feature obviously cannot be simply due to 
a possible miscalculation (underestimate) of 6, the back- 
ground correction. Indeed, by underestimating b, we could 
find false emission-like features at wavelength where the 
flux is low. But, one should increase the estimate of b 
above its 4tT upper limit (8 • 10~^^ erg cm~^ s""'^ A~^) to 
explain the detected emission feature only by this effect. 
In addition, the emission is detected in both Fe ii lines, 
in particular in the Fe ii line at 2374.461 A where the 
level is far above the zero level, and for which an error 
on the estimate of b has almost no effect on the result. 
In fact, b is well-determined by the very bottom of the 
strongest Fe ii lines where the level is clearly less than 
1 • IQ-" erg cm-2 s'^ k'^ (Fig. |l|). 

However, apparent emissions due to the addition of 
the statistical noise and a bad estimate of the background 
level might be possible. But, although not excluded, it is 
very unlikely that this coincidence can contribute to give 
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Wavelength (A) 

Fig. 2. Ratio of the disk spectrum to the star spectrum 
(top histogram, left axis) and plot of the star spectrum #2 
(bottom spectrum, right axis). The pixels have been re- 
binned in boxes of 50 kms~^width. To evaluate the ratio, 
the disk spectrum has been rescaled with = 197 and 
6 = 3- 10^^'^ erg cm^^ s^^ A^^. The crrorbars are at la- 
level. The points from which a signal is detected at 2a are 
plotted in bold. The dip in the spectrum between 2377 A 
and 2378 A is due to hot diodes. An excess emission is 
clearly detected in the red part of the two strongest Fe ii 
lines. 



apparent emissions in the two sides of the two strongest 
Fe II lines. 

The most important alternative to emission by Fe ii is 
a time variation of the /3Pictoris spectrum between the ob- 
servations of the template (the star) and the disk spectra. 
If the absorption component in the Fe ii lines had signif- 
icantly decreased during the acquisition of the data, then 
the result is an apparent excess of emission in the second 
spectrum obtained with the slit off the star. Although the 
time between both spectra has been minimized, this pos- 
sibility cannot be excluded without new observations, for 
example on the other side of the disk where the emission 
should be stronger in the blue. 
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Fig. 3. Plot of the Fe ii emission in the disk spectrum 
at 0.5" obtained by the subtraction of the star scattered 
light (top histogram, solid line for the 2374.461 A Fe ii 
line and dashed line for the 2382.765 A Fe ii line, left 
axis). The data have been rebinned in boxes of 50 kms~^ 
wide. Abscissa are radial velocities relative to the radial 
velocity of the star, i.e. heliocentric minus 21 kms~^. The 
errorbars are at la level. The points for which a signal 
is detected at 2a are plotted in bold. The two bottom 
spectra are (3 Pictoris spectra (right axis). Because of the 
strong Fe ii absorption lines of the stable disk (at kms^^ 
relative to /3 Pictoris, 21 kms~^ heliocentric), pixels with 
a velocity between —25 kms^^ and -f 20 kms^^ relative to 
(3 Pictoris are very noisy and are not taken into account. 
The two Fe ii lines give about the same emission although 
their shapes and their levels are very different. 



3. The Fe ii emission lines 

3.1. The dynamics of the Fe ii ions 

If this detection is really due to emission by Fe ii ions, the 
lines width must be explained through the dynamics of 
these Fe ii ions in the disk. The Fe ii ions must be ejected 
from the /3 Pictoris system by the radiation pressure which 
is stronger than the gravitation by a factor PfcII ~ 5 
(Lagrange et al., 1998). After ejection, they rapidly reach 
a constant asymptotic velocity Uoo- If they are ejected from 
a body on a circular orbit, Voo ^ \J {2(3 — l)(GM/ao), 
where qq is the radius of the orbit of the parent body. 
If they are ejected from a comet on a parabolic orbit, 
the final velocity is Uqo ^ \/ {2(3GM / ao). In this simple 
scheme, the observed final velocity of about 100 km s~^ 
(Fig. H) corresponds to a production at about 1.5 AU from 



the star, or similarly to the absence of gas drag beyond 
that distance. 

The emission lines are stronger in the red than in the 
blue. This is similar to the asymmetry already observed in 
the cometary absorption lines which are mainly redshifted 
(Beust et al. 1996). This last asymmetry is well- explained 
by the evaporation of comets with a small range of lon- 
gitude of periastron (Beust et al. 1998). An alternative 
explanation for the observed asymmetry in the emission 
lines could be an extended shape for the dragging torus 
of gas needed to support the radiation pressure on the 
Ca II and Fe ii ions observed at zero radial velocity in 
the stable gaseous disk (Lagrange et al. 1998, Beust et al. 
1998). For both explanations, it is clear that the observed 
south-west branch of the disk must be the "red side" of 
the disk (Fig. |). This provides an observational test : the 
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north-cast branch must present a larger emission in the 
blue lines. 



3.2. The Fe ii density 

The total brightness of Fe ii emission lines can be evalu- 
ated to be: 



emission 



4tt 



n{r)/r^dx 



(3) 



in erg cm ^ s ^, where is the solid angle covered by 
the spectrograph slit. The SSA slit (0.22"x0.22") gives 
fl = 10""*^^. d is the distance to /? Pictoris {d = 19.3 pc« 
6 • lO^^cm). s is the frequency integral of the cross section. 
F^^^'^ is the brightness per unit of frequency of /3 Pictoris 
seen from the Earth at the relevant wavelength. n(r) is 
the density of the observed ion at a distance r from the 
central star, dx is the differential length along the line of 
sight. We can define a weighted integral equivalent to the 
column density by 



N, 



ro 



_ f n{r) 2 



47rr2 F ■ ■ 
7 ' emission 

Qd^sFr'' ' 



(4) 



where ro is the impact parameter of the line of sight 
(ro(0.5") = 10 AU). For the Fe ii 2382 A line, s = 
8-10- 

Ff ^'^AA ~ 10" 



^ cm^ s ^ and the observed emission is F^rmssion 
^^erg cm~^ s~^. Finally, we have 



A^ioAu = 5.2 • 10' 



,26 



^ emission 

erg cm~2 s~^ 



cm 



-/VioAu ~ 5 • lO^^cm 



(5) 



(6) 



This value is consistent with the Fe ii column density 
(^Fe II ~ J n{r)dr — 3 • 10^** cm^^) and the hypothesis 
that Fe ii should be gathered in the dragging torus around 
1 AU and that they have an r^^ distribution beyond this 
torus (Lagrange et al. 1998). To evaluate the Fe ii volume 
density, no, at ro, we can define the dimensionless quantity 
K hy K = J 'n{r)rQ / {nor^)dx. Then, 



no 



AirroF, 



emission 



N 



nd'sFr'^K roK 



(7) 



If we make the assumption that n(r) = no(ro/r)", we 
have Ka=o = tt, Ka=i = 2, Ka=2 = Tr/2, K^^s = 
357r/128, et cetera, a = 2 would correspond to gas expelled 
by radiation pressure. For reasonable value of a, within a 
factor of 2, ii' ~ 2. As a final result, we get 



no ~ 2 ■ 10 



at ~ 10 AU. 



(8) 



The calculation has been done with the hypothesis of 

an optically thin line. We also assume that the filling factor 
in the vertical direction is 1, thus the obtained density no 
is in fact a lower limit. The other Fe ii line at 2374 A has 
an oscillator strength 10 times smaller but the stellar flux 
is larger at this wavelength (~ 4-10""'^^ erg cm~^ s~^ A~^). 
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Scattering 
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Fig. 4. Sketch of the possible dynamics of the P Pictoris 
disk. Fell ions are produced close to the star, spiral into 
a massive torus at about 1 AU and then cross the line 
of sight defined by the slit. We see the resonant scatter- 
ing through emission lines blueshifted and redshifted by 
the radiation pressure. The observed larger amount of red- 
shifted Fe II is due either to the rotation of the massive 
torus which drags these Fe li ions, or to an asymmetri- 
cal production of ions as inferred from the falling comets 
scenario. In summary, the large velocities are due to the 
radiation pressure, and the asymmetry between the red 
and the blue is due to the rotation of the system. 

It gives about the same value of density within a factor 
of 2. 

The corresponding Fe ii production rate can be 
roughly estimated by the flux of material (v — 100 km s~^) 
through the area observed at 10 AU with high ^"(0.22") = 
4.4 AU: Q(Fe ii) = 2'KroHvno « 10^ kg s^^. With solar 
abundances, this corresponds to the total disruption of 
about 10 asteroids per year with a radius of 30 km. 

4. Need for confirmation 

To confirm the tentative detection presented in this pa- 
per, new observations are really needed. The first obvi- 
ous method will be to observe the other (north-east) side 
of the disk which should present the same feature ex- 
cept blueshifted instead of redshifted. We can also observe 
other lines which should present the same characteristics. 
The Fe ii lines at 2600 A and Mg il lines at 2800 A arc 
well-suited candidates. The Cr ii line at 2050 A with an 
intermediate jj « 3 is also an interesting target (La- 
grange et al. 1998). STIS, with long-slit capability is tech- 
nically better suited to this observation. 
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5. Conclusion 

HST observations planned to detect the emission lines of 
ions in the /3Pictoris disk with the spectrograph slit placed 
off the star gave a marginal detection of possible emission 
lines from Fe ii ions at 0.5" from /3 Pictoris. This would 
suggest an Fe ii density of ~ 2 • 10~^cm~^ at 10 AU. 

If real, this indicates a large production rate of gas 
and dust equivalent to the total disruption of 10 bodies 
of 30 kilometers in radius per year. This corresponds to a 
production rate of about 2 • 10~''MEarth per year. New ob- 
servations are obviously needed to confirm this detection. 
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